Spin-spiral states in undoped manganites 
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The experimental observation of multiferroic behavior in perovskite manganites with a spiral spin 
structure demands to clarify the origin of these magnetic states and their relation to ferroelectricity. 
We show that spin-spiral phases with diagonal wavevector and also the E-type phase exist for 
intermediate values of the Hund's rule and the Jahn- Teller coupling in the orbitally ordered and 
insulating state of the standard two-band model Hamiltonian for manganites. Our results support 
the spin-current mechanism for ferroelectricity and present an alternative view to earlier conclusions 
where frustrating superexchange couplings were crucial to obtain spin-spiral states. 
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The coexistence of long-range magnetic order with 
spontaneous electric polarization is commonly refered to 
as multiferroic behavior In recent years multiferroic 
materials have attracted special attention from the con- 
densed matter community because of their potential for 
applications in memory and data storage devices @, Q ■ 
Despite the initial observation that materials with coex- 
isting ferroelectric and magnetic orders are rare in nature 
[H, an increasing number of multiferroic materials with 
interesting properties has been discovered Differ- 
ent mechanisms for the origin of multiferroic behavior 
have been proposed and partially identified [9j-fl4j. In a 
selected class of multiferroics ferroelectricity is driven by 
the existence of a non-trivial magnetic order, e.g., with a 
spiral spin structure [l5|-[l8j]. It has therefore become a 
key issue to clarify the conditions under which non-trivial 
spin states can occur in different models and materials. 

Hole-doped perovskite manganites are known for their 
rich phase diagrams and complex transport phenom- 
ena [19( | . The recent observations of ferroelectricity in 
TbMnC>3 and DyMnC>3 stimulated further research also 
on the undoped materials. The magnetic groundstatc 
of RM11O3 with R = La, Pr, Nd and Sm is an A-typc 
antiferromagnet (AFM), with ferromagnetic (FM) order 
in the a-b plane and a staggered spin pattern along the 
c-axis. It changes to a spiral magnet for R = Tb and 
Dy, and finally to an E-type AFM for R = Ho, where 
zigzag FM chains alternate in their preferred spin direc- 
tion [T7l . [20j . Although the A- type order of the proto- 
type compound LaMnC>3 is well understood in terms of 
Goodenough-Kanamori rules and orbital ordering fl9j |. 
the magnetic structure of the materials with smaller ionic 
radii are much less understood. It was proposed that due 
to GdFe03-type structural distortions longer range inter- 
actions become relevant and thereby may lead to com- 



plex magnet ic g roundstates with spiral or E-type spin 
patterns [13, Ef. The existence of an E-type pattern in 
these models with longer range interactions has been put 
into question recently p2j . 

In this letter we explore the magnetic groundstates 
in a two-band model for RM11O3 without invoking the 
next-nearest neighbor or even longer range interactions. 
The model consists of itinerant electrons coupled locally 
to core spins via the Hund's rule coupling Jh, and a 
nearest neighbor antifcrromagnetic interaction between 
the spins. Importantly, we refrain from taking the com- 
monly adopted double- exchange limit Jh — > 00. The 
two-orbital nature of the model allows for the inclusion 
of the Jahn- Teller (JT) distortions as a source for or- 
bital order and hence the insulating character of the un- 
doped manganites. We find that both the spiral spin 
states and the E-type states are stable in the parameter 
regime relevant for manganites. The experimentally ob- 
served transitions between different magnetic states are 
obtained only for finite JT couplings and intermediate 
values of Jh- 

Specifically, we consider a two-dimensional two-band 
model with the Hamiltonian 1231: 
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at quarter filling to describe the undoped manganites 
RMJ1O3. Here, c iaa and cJ Qcr are annihilation and crc- 

or- 



ation operators for electrons with spin a in the e g 
bital a € {x 2 — y 2 ,3z 2 — r 2 }, which from here onward; 
is labeled as 1 and 2, respectively, t^g denote the hop* 



2 



ping amplitudes between the two e g orbitals on nearest- 
neighbor sites and have the cubic perovskite specific form 
[if: tf, = t y n = t, q 2 = t\ 2 = t/3, tf 2 = t* 21 = -t/y/3, 



t v 



t 21 = t/^/3, where x and y mark the spatial direc- 



tions. In a commonly used approximation for mangan- 
ites, the core spins are treated as classical vectors with 
|S| = 1; the justification of this approximation was quan- 
titatively verified [24[ ■ tr ia denotes the electronic spin op- 



erator defined as er^ 



,t 



a' C iaa L ao' Cia °' 



where r M are 



the Pauli matrices. The one-band version of the above 
model and the Kondo lattice model with J s = have 
been previously analyzed in one- and two-dimensions in 
search for non-trivial magnetic groundstates (25l - [28j . Al- 
though spin-spiral states were found for selected com- 
binations of band fillings and Hund's rule coupling, the 
connection to the spin spirals observed in the RMn03 is 
unclear since the insulating character of orbitally ordered 
RMn03 is not captured by the one-band models. 

By applying a canonical transformation we rewrite the 
Hamiltonian in a basis where the spin-quantization axis 
is site-dependent and points along the direction of the 
local core spin. Introducing polar and azimuthal angles 
9 and 0, respectively, the transformation is defined as: 
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FIG. 1: The cone angle 6 and the wavevector (q x , q y ) of the 
lowest-energy spiral state as a function of Jh for, (a) J s — 
and (b) J s = 0.04£. The inset in (a) shows the J s dependence 
of the critical value of Jh required to obtain the FM state. 



For this choice of variational spin states the Hamilto 
nian matrix reduces to a 4 x 4 matrix after Fourier trans 
formation with, d las = N' 1 / 2 ^ k e -ik 
kinetic energy term in Eq. (3) is then written as, Hk 

J2 k r>l H(k) D k , where, D* = ~ 

and the 4x4 matrix % is given by 
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Here d iap (d iaa ) annihilates an electron at site i in orbital 
a with spin parallel (antiparallel) to the core spin. In 
terms of d operators the Hamiltonian reads: 

(ij)a s,s' 

^ ^ { n iap n iaa) Js ^ ] Sj • Sj , (3) 

ia (ij) 



with n ias 
given by, 



d ias d ias . The coefficients f ss t, are explicitly 
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The advantage of this transformation is that the Hund's 
rule term now becomes diagonal. For a fixed configura- 
tion of classical spins the Hamiltonian is bilinear in the 
fermion operators. However the one-particle Schrbdinger 
equation can not be solved in closed form for an arbitrary 
spin configuration. We therefore selectively analyze the 
core spin configurations described by the polar and az- 
imuthal angles 9i = 6 and <f>i = q ■ . We refer to as 
the cone angle and q as the spiral wavevector. Despite 
this restriction most of the well known spin patterns as 
observed in various magnetic materials are included, e.g., 
ferromagnetic, antifcrromagnctic, cantcd-fcrromagnctic, 
and spin-spiral patterns. 
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with the matrix elements, 
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with = k^ ± Qfi/2. For a single band Hkin reduces to 
a 2 x 2 matrix structure and therefore the eigenspectrum 
is straightforwardly obtained in a closed form expression 
[2l| . For the two-band case considered here the closed 
form result for the dispersion is rather involved and we 
therefore diagonalize the above 4x4 matrix numerically 
for each momentum k on finite lattices of upto 256 x 256 
sites. 

In Fig. 1 we plot the values of the cone angle O and 
spiral wavevector (q x , q y ) corresponding to the minimum- 
energy state as a function of Jh- Over almost entire 
parameter regime the cone angle is found to be 7r/2, 
which corresponds to planer spin states. Another fea- 
ture is that q x = q y , suggesting that the diagonal spi- 
rals are more stable, which is consistent with the ex- 
perimental findings [20|. For the diagonal spirals with 
(ixiQy) = the spiral pitch q smoothly vanishes 

upon reaching the FM state with q = 0. Close to the FM 
phase the cone angle of the spiral state slightly deviates 
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from 7r/2. The grovmdstate jumps discontinuously from 
an antifcrromagnct with q — tt to a spiral with q < it. 
Inset in Fig. 1 (a) shows the variation with J s of the crit- 
ical value of Hund's rule coupling J H for the transition 
to a ferromagnet. 

Before arriving at the magnetic phase diagram of the 
Hamiltonian, it is essential to consider other states which 
arc not captured by the ansatz (#,, </>i)=(6, q-Tj) [26l.[29j. 
The energies of the earlier suggested candidate states, in- 
cluding the E-type states, for the magnetic order in man- 
ganites are obtained by exact numerical diagonalization 
and compared with those of the spiral states. 

The groundstatc-phasc diagram for the quarter filled 
system is shown in Fig. 2(a). The E-typc phase is stable 
in a wide region of parameter space. Spiral states are 
favored for larger values of J s , and also in a narrow win- 
dow between the FM and the E-type states. Here, the 
FM state in a two-dimensional model is representative of 
a single plane of the A-type AFM state. In the experi- 
ments on RMn03 two transitions arc observed, first from 
the A-type AFM to the spiral state and subsequently to 
the E-type phase, upon reducing the ionic radius of the 
rare-earth element [17|- Therefore the existence of a di- 
rect FM to E-type transition appears to contradict the 
experimental observation. Moreover, the phase diagram 
of Fig. 2(a) corresponds to a metallic state with a finite 
density of states at the Fermi level, the undoped man- 
ganites however are insulating. Therefore it is essential 
to include the source for the opening of an energy gap 
in the spectrum in order to obtain results applicable to 
RMn03. For this purpose we add an adiabatic Jahn- 
Teller coupling to the Hamiltonian, which is given by, 
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In Eq. (7), Q X i and Q z i are lattice distortions corre- 
sponding to two different JT modes. t x % = X^o-( c !ic 



4<7 C ila) and T zi = XL 
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i2a C i2a) are Orbital 



pseudospin operators [19|. In the undoped manganitcs 
staggered JT distortions are accompanied by orbital or- 
dering with transition tempertatures much hig her than 
the temperatures scale for magnetic ordering [17| . There- 
fore the pattern for the JT distortions is expected to be 
robust upon cooling even though the magnitude of the 
distortions may depend on the magnetic structure. A 
real-space Monte Carlo study verified the staggered or- 
dering of the Q x component when Hjt is included in 
the Hamiltonian [30l |. Therefore we adopt this pattern 
for the lattice distortions and set Q X i = Q x e 1 v ,r > 7r )' r< and 
Q z i = 0. The second term in Eq. (7) is the energy 
cost associated with the distortion of the lattice with 
|Qi| 2 = Q xi + Q 2 zi . We set the stiffness constant K = t 
as in previous theoretical model analyses of manganitcs 
[l9l . Ho] • We treat Q x as a variational parameter in the 
calculations and optimize it by minimizing the total en- 
ergy. Therefore Q x is allowed to vary in the different 
magnetic states. 
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FIG. 2: (Color online) Groundstate-phase diagrams J a ver- 
sus Jh in the absence (a) and presence (b) of Jahn- Teller 
distortions. The spiral states have the cone angle = n/2 
except in the narrow region between FM and the E-type states 
in (a). In both cases q x — q y reduces monotonically as the 
FM state is approached. 



Due to the staggered orbital order the determination 
of the eigenspectrum becomes slightly more involved. 
Instead of a 4 x 4 matrix, we now have to diagonal- 
~ HQs) M 
M -H(k') 



ize the 8x8 matrix 
k' = k + (it, 7r) and 



for each k with 



M = 
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The energy minimization procedure with respect to 
and q is followed as before but in the presence of stag- 
gered Jahn- Teller distortions, which lead to an insulating 
state with staggered orbital order. The resulting phase 
diagram is shown in Fig. 2(b). Remarkably the inclusion 
of the Jahn- Teller distortions leads to the appearance of 
spiral states in between the E-type and the FM phases 
in a wide range of the Hund's rule coupling. 

An important effect of the size reduction of the rare 
earth ion in RMn03 is the decrease of the electronic 
bandwidth due to the Mn-O-Mn bond angle moving fur- 
ther away from 180°. For the model calculations, it is 
simpler to vary J s and A rather than changing the hop- 
ping parameters which control the bandwidth. In Fig. 3 
we therefore show the phase diagrams in the parameter 
space of J s and A for two representative values of Jh- 
For Jh = 5t, the small J s regime is ferromagnetic, which 
describes a single plane of the A-type AFM observed in 
RMnOs with R = La, Pr, Nd and Sm. A two-step tran- 
sition occurs from FM via the spiral to the E-type state 
by increasing both A and J s (indicated by the arrow in 
Fig. 3(a)), which effectively translates to reducing the 
bandwidth. In the (1,1) spiral state the pitch q increases 
along the direction of the arrow (see Fig. 3(a)). The val- 
ues of q at the two end-points of the planar spiral phase 
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FIG. 3: (Color online) J s -A phase diagrams for (a) Jh = 5t, 
and (b) Jh = lOOi. The arrow in panel (a) is indicative of the 
path traced upon reducing the ionic radius of the rare-earth 
ion in RMn03, which leads to the correct sequence of tran- 
sitions from FM to spiral to E-type states. Such a sequence 
can not be traced along any straight line in panel (b). 

along the arrow are 0.137T and 0.25-7T. These values are 
close to the experimental results of 0.147T and 0.197T ob- 
tained for TbMn0 3 and DyMn0 3 , respectively 0. The 
pitch vector q' = (0, fc,0) reported in the experiments 
translates to q = (fc/2,fc/2,0) on the Mn square lattice 
used here due to a 45° rotation between the two coordi- 
nate systems. The E-type AFM state has been observed 
in HoMn0 3 [13, ■ 

For the larger Hund's rule coupling Jh the stability 
region of the spiral states shrinks considerably (see Fig. 
3(b)). Moreover, the spiral states are no longer sand- 
wiched between the FM and the E-type phases. This 
implies that it is not possible in the commonly adopted 
double-exchange (Jh — > oo) limit to find transitions to 
spiral states as observed in RMn0 3 upon varying the 



bandwidth, unless further interactions are added to the 
two-band model Hamiltonian. 

The presence of a spiral structure in an insulator has 
been identified as one possible source for a spontaneous 
electric polarization P by generating spin currents [Til ] . 
The direction of P is perpendicular to both the direction 
of the spiral pitch vector q and the cone axis of the spi- 
ral (l2| . Our model of choice is isotropic in spin space 
and thus can not determine the orientation of the cone 
axis relative to the crystallographic directions. Using the 
input from the experiments regarding the direction of 
the cone axis for the spiral state we indeed obtain the 
direction of P consistent with the experiments 0, [3l| . 
Within the spin-current mechanism the magnitude of P 
is controlled by the length of the pitch vector q, which 
we obtained in the experimentally relevant range. 

Therefore our results verify that already the standard 
two-band model for manganites has all the necessary in- 
gredients to sustain the magnetic spiral and the E-type 
phases as observed in the undoped perovskite mangan- 
ites. A finite Hund's rule coupling of the order of the 
bandwidth leads to a spiral pattern and a wavelength 
which both compare well with the observed magnetic 
structure in TbMn0 3 and DyMn0 3 . These results sup- 
port the applicability of the spin current mechanism as 
the source for ferroelectricity in RMn0 3 . Lattice distor- 
tions of the GdFe0 3 -type are likely to give rise to addi- 
tional longer-range couplings, which may further stabilize 
the spiral and the E-type states. The essential physics of 
these non-trivial spin states is already contained in the 
simpler short-range two-orbital model with finite Jh- 
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